Transport through the single- molecular dots in an external irradiation 
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We present a fully nonequilibrium calculation of the low-temperature transport properties of a 
single molecular quantum dot coupled to local phonon mode when an ac field is applied to the gate. 
The resonant behavior is shown in the time-averaged differential conductance as the ac frequency 
matches the frequency of the local phonon mode, which is a direct consequence of the satellite- 
phonon-peak structure in the dot electron spectral function. The different step structure with and 
without the external irradiation is found in the I-V curves, and the oscillation behavior is found in 
the step height as a function of the irradiation intensity. 

PACS numbers: 73.40.Gk, 73.63. Kv, 85.65.+h 
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Recent advances in nanotechnology have allowed the 
fabrication of very small molecular quantum dots weakly 
coupled to the macroscopic charge reservoirs (leads) 0. 
In contrast to the semiconductor dot, which is quite rigid 
in space, the molecules involved in the electron tunnel- 
ing process naturally possess the vibrational degrees of 
freedom which will inevitably react to the transform of 
electrons through the molecular quantum dotsQ. In 
addition to the importance in molecular-scale electron- 
ics from the application point of view, these artificial, 
tunable devices are potentically important for under- 
standing the basic physics including the many-body ef- 
fect. Theoretically, a lot of effort has been focused on 
the quantum conductance of molecular systems based 
on the kinetic equation approach|2j, the rate equation 
approach^, the correlation effects |J, the nonequilibrium 
quantum theory [fj |f| 0, El, and the strong coupling to 
environment 0. So far, the stationary quantum trans- 
port through the molecular dots has been considered, 
while the influence of a time-dependent ac field on the 
current has not been well addressed. Irradiation of a 
quantum dot with an ac offers a new way of 

affecting its dynamics, which enables one to study the ef- 
fect of electron-phonon interaction on the transport phe- 
nomenon of molecular dots in essentially nonequilibrium 
condition. 

In this Letter, we use the Keldysh nonequilibrium 
Green's function technique to study the nonlinear ac 
transport through a single-molecular quantum dot cou- 
pled to a local phonon mode with the external irradi- 
ation applied to the gate, for the first time. After a 
canonical transformation, we obtain a formula for the 
time-dependent current in general terms of bias, tem- 
perature, the intensity and frequency of the external ac 
field, and the electron-phonon coupling. We show that 
the satellite-peak structure due to the electron-phonon 
interaction can be probed by imposing on top of gate 
bias an ac bias voltage. The satellite-peak structure in 
the dot electron spectral function gives rise to resonant 
behavior in the time-averaged current as the ac frequency 
matches the frequency of local phonon mode, which can 



be observed directly in experiments. The calculated I- 
V curves also show the different step structures with 
and without the external irradiation, and the step height 
shows the Bessel-type oscillation behavior as a function 
of irradiation intensity. 

In this work we consider a simplest Holstein-type 
model with a single phonon mode is employed to ad- 
dress the vibrational degrees of freedom in the molecular 
dot. All other complexity of real molecular devices, apart 
from interaction with phonon mode, is ignored. Then the 
system Hamiltonian can be written as 



H = H, 



leads 
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c L;cr ( c kri<j) arc creation (annihilation) operators for the 
noninteracting electrons with momentum k and spin in- 
dex a in the left (77 = L) or right (77 = R) metallic leads. 
<x>o is the frequency of the single phonon mode, and a' (a) 
is the phonon creation (annihilation) operator. Ho de- 
scribes the electron in the quantum dot coupled to the lo- 
cal phonon mode with the coupling constant A, where d' a 
(d a ) is the dot-electron creation (annihilation) operator, 
and eo (t) is the single energy level of the dot which can be 
tuned by the external irradiation, eo (t) = eo+V r cos (ui r t) 
for harmonic bias. Here we assume that the metallic leads 
are dc biased, neglecting the possible "leakage" of the ir- 
radiating ac field to the leads. The generalization onto 
the case of nonzero ac bias is straightforward. Ht de- 
scribes the tunneling coupling between the dot and the 
leads, where the tunneling matrix elements V£ a a transfer 
electrons through an insulating barrier out of the dot. 
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Based on the Keldysh noncquilibrium Green's func- 
tion formarism|l2j , the time-dependent current from the 
X] lead to the dot is given bv[lO| 

vo- -f {*■•<•-•'>., 

x [G< a {t,t') + .fn^)G r aa {t,t')]}, (3) 

where Jl(b) ( e ) arc the Fermi distribution function of the 
left (right) leads, which have different chemical poten- 
tials upon a dc bias voltage /jll — )JtR = eV . T 7 ^ = 
2 

2-Kpa (0) VP a a characterizes the coupling between the 

dot and the leads, and p a (0) is the spin-a band density 
of states in the leads. Here we assume that the leads give 
rise to a flat, energy independent, density of states (i. 
c, the wide-band limit). G rty< ^ is the retarded (lesser) 
Green's function of the dot. 

In order to compute the time-dependent current, one 
has to compute the dot electron Green's functions in the 
presence of both the electron-phonon interaction and the 
tunneling coupling between dot and leads. The Green's 
function can be calculated by performing the canonical 
transformation S = (X/ujq) ^2 a dt a d a (at — ».) (lflj. and 
then the dot level is renormalized to eo — A, where 
A = X 2 /u>o, & n d the tunneling coupling term is also 



VP c 



d a X + H.C. 



renormalized as Ht = J2k n a a 

where X = exp [— (X/wq) (at — a)] . Ignoring the 
effects of narrowing the bands of leads due to the 
phononsflif. the dot-electron Green's function can be 
decoupled as G r aa , (t,f) = G r aa , (t, f) (X (t) (t')) ph , 

where & aa , (t,t>) = -iQ (t - it) ({rf a (t) , dt a (t>)}) ^ 

d a (t) 
X(t) 
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H e t = Hx, and 
e iH ph tx e -%H ph t_ renorma ii za tion 
factor due to the electron-phonon interaction 
isd (X(t)XHt')) = e -*( 4 -0, where <t> (t) = 



(A/w ) [Nph (1 - e 
N ph = 1/ [exp(/3w ) 



' ph 



) + (N ph + l)(l-e-' l ^% and 



I]. The retarded Green's function 
b; 

and the result is 



can be easily obtained by the standard Dyson equation 
approach yjj, 



Gw(M') - -i6 aa ,0(t-t')e- l ti dTV " cos ^T) 

Xe -i[ eo -A-jr a ](i-t')-*(i-t') j ^ 

where T a = r^+T^ is the total tunneling coupling to the 
leads. Following the operational rules |l5| to the Dyson 
equation for the contour-ordered Green's function, the 
Keldysh Green's function is found to be 



(*,*') = <W / dh j dt 2 G r aa (tM) 

xY,<(h,t 2 )G a aa {t 2 ,it), (5) 



with the less self-energy 

S< (*!,**)=*£ I grSe-»[-^-*»)]/ fI (e)e- fa (*-*-). 
i 

(6) 

Without electron-phonon interaction, the above result 
fully agrees with that for time-dependent transport 
through nonintcracting quantum dot|lfl|. 

Substitution of the Green's functions into Eq. gives 
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x / dt 1 e- T ^ t - tl '>Jm[A a (u,t l )]\, (7) 



mthA a {u,t) = J^dt'e^-^G^it^'). Obviously, 
in the time-independent case, A a (to) is just the Fourier 
transform of the retarded Green's function G r aa (10). Af- 
ter some algebra, we find that for this model, 
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i(l — m)u} r tnuJo{3/2 



[2g[N ph (N ph + l)] 1/2 } 



lu — (eo — A) + moj r — tiluq + ^T a 



(8) 



where J m (z) is the Bessel function of the m-th order, 
/„ (z) is the Bessel function of complex argument, and 
I, to, n = 0, ±1, ±2---. Eq. © together with the 
current expression Eq. (JJJ provides the complete solu- 
tion to the time-dependent transport of molecular quan- 
tum dot coupled to local phonon mode in the ac field. 
Experimentally what is interest is the current on a time 
scale long compared to 2tt /tu r . Here we discuss the time- 
averaged current (J (t)), which could be directly relevant 
to experiment. For this model, we then obtain 

(J L (t)) = -(J R (t)) 
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(9) 



In the absence of external irradiation, i.e., V r = 0, Eq. 
© fully agrees with the result of dc current in Rcfs. 
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FIG. 1: The time-averaged differential conductance, in units 
of 2e 2 /h through the molecular quantum dot as a function of 
the Fermi energy Ef measured relative to the single level of 
the dot in the absence (V r = 0), and in the presence of ac 
field with the frequency ui r = <^o and three different values 
of irradiation intensities V r = 0.5u>o, ^>o, and 1.5ojo, respec- 
tively. The energy is measured in units of the frequency of 
the phonon mode u>o, X = OAujq, eo = 0, and T — 0.2ljo. 
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FIG. 2: The time-averaged differential conductance, in units 
of 2e 2 /h as a function of the frequency of irradiating ac field 
for different intensities of ac field with V r /u>o = 0.25, 0.5, 1, 
and 2, respectively. Here X/coo = 1, eo = 0, T = 0.2wo, and 
Ef = (-l)wo- 



P, |8| . The contribution of the external irradiation to the 
transport becomes significant when the argument of the 
Bessel function, V r /uj r , is of the order of unity. Another 
important consequence is that the time-averaged current 
under irradiation is proportional to J^. As a result, one 
should expect the change in current due to the external 
irradiation to depend on the intensity of irradiation. 



FIG. 3: Current- voltage characteristics, in units of 2e/h, 
without (V r = 0) and with the ac field for the different ir- 
radiative frequency: uo r = 0.5cjo and lwo- Here, X/luq — 1, 
Ef/ujo = -1, e = 0, Vr/coo = 2, and T = 0.05w . 



For simplicity, we consider the tunneling coupling be- 
tween the molecular dot and the two leads to be symmet- 
ric and independent of the spin index, i.e., = = 
= = r. In Fig. 1, we plot the zero-temperature 
differential conductance as a function of the Fermi en- 
ergy measured relative to the single level eo of the dot 
in the external ac field with frequency uj r = uq and dif- 
ferent values of irradiation intensities. For comparison, 
we also plot the differential conductance in the absence 
of ac field. At zero ac bias voltage, our results agrees 
well with that in Refs. 0, 13 > where the electron-phonon 
coupling can lead to the satellite resonant peaks. Fig. 
1 shows that the ac field with frequency Lu r = ujq can 
lead to the enhancement of satellite resonant peaks in 
the positive energy region and the appearance of new 
peaks in the negative energy region. One can also see 
that the main peak is suppressed by the irradiation while 
increasing the intensity. In Fig. 2, we plot the zero- 
temperature differential conductance as a function of the 
frequency of ac field for different values of irradiation in- 
tensities by fixing the Fermi energy of the leads Ep as 
Ef — (eo ~ A) = to avoid unnecessary complications. 
As the intensity of ac field increases, resonant signals 
are clearly shown when the frequency of ac field satisfies 
muj r = tiluq. The satcllitc-phonon-pcak structure in the 
dot electron spectral function, as shown in Refs. 0, Q 
and also in Fig. 1 of this Letter, gives rise to resonant be- 
havior in the conductance as the irradiating ac frequency 
matches the frequency of local phonon mode, and can be 
observed directly in experiments. Fig. 2 also shows that 
more resonant signals can be observed while increasing 
the irradiation intensity. 

Fig. 3 shows the calculated zero-temperature current- 
voltage curves with and without the external ac irradi- 
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FIG. 4: The height of the zeroth (N = 0) [a] and the first 
(TV = 1) [b] step, in units of 2e/h, as a function of the irradi- 
ation intensity by fixing the Fermi energy of the leads Ef as 
Ef — (eo — A) = 0. Here, ui r — wo, and T = 0.05cjq. 



ation. Here we assume the leads to be symmetrically 
voltage biased, i.e., +V/2 on the left lead and — V/2 on 
the right one, to avoid unnecessary complications. In the 
absence of external ac field, clear steps appear at roughly 
2wo intervals in the weak tunneling coupling limit, cor- 
responding to eV/2 = ±nw ; with n = 0, 1, 2 ■ • • , and 
the height of the N-th step decreases with N, which can 
be easily understood from Eq. (|10f> [see below] by taking 
V r = for small electron-phonon interaction. In the pres- 
ence of ac field with frequency co r = too, steps appear at 
the same intervals 2loq as those in the zero ac field, while 
the step height is modulated by Bessel function due to 
the irradiation [see Eq. (fLU|) in below] . Fig. 3 also shows 
that more steps appear at the intervals loq in the case of 
u> r = 0.5lu$, corresponding to eV/2 = muj r — ntuo with 
m = 0, ±1, ±2 • • • . The zero-temperature N-th step's 
height in an external ac field with frequency to r = loq can 
be analytically obtained from Eq. @ at the fixed Fermi 



energy Ep as Ep — (e — A) = 0, 
AJjv 



/2e\ r^rg 



n=0 



T 2 

J n+N 



/2 



N 



V r 



(10) 



In Fig. 4, we plot the height of the N-th. step, where 
N = and 1, as a function of the irradiation inten- 
sity for different values of the electron-phonon coupling 
constant. The oscillation behavior of the step height is 
clearly observed for the small electron-phonon interaction 
due to the external irradiation [see Eq. (|10fl ], while this 
oscillation smears out for large electron-phonon interac- 
tion. The large electron-phonon interaction enhances the 
processes of absorption and emission of many (n > 1) 
phonons when electron tunnels through the dot, and then 
the summation of Bessel function with large values of in- 
dices n ± N results in the smearing of the oscillation 
behavior [see Eq. 1)10(1 ], 

Summarizing, using the Keldysh nonequilibrium 
Green's function technique, we have studied in this work, 
to the best of our knowledge, for the first time the time- 
dependent transport through a single molecular quan- 
tum dot coupled to a local phonon mode in the presence 
of an external ac field. We show that the external ir- 
radiation provides another important experimental tool 
where both the equilibrium and out of equilibrium trans- 
port phenomenon can be probed. In particular, reso- 
nant behavior as the ac frequency matches the frequency 
of local phonon mode is shown to exist as a result of 
the satellite-phonon-peak structure in the dot electron 
density of states. The nonlinear I-V curves exhibit new 
structure caused by the external irradiation, which can 
be investigated experimentally. 
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lead to this work and Z.-R. Liu for useful discussions. 
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